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Abstract 



Kumar et al. (2006) obtained a fifth order polynomial in lu for the dispersion relation 
and pointed out that the calculations preformed by Porter et al. (1994) and by Dwivedi & 
Pandey (2003) seem to be in error, as they obtained a sixth order polynomial. The energy 
equation of Dwivedi & Pandey (2003) was dimensionally wrong. Dwivedi & Pandey (2006) 
corrected the energy equation and still claimed that the dispersion relation must be a sixth 
order polynomial. The equations (11) — (19) of Dwivedi & Pandey (2006) and the equations 
(24) — (32) Kumar et al. (2006) are the same. This fact has been expressed by Kumar 
et al. (2006) themselves. Even then they tried to show this set of equations on one side 
gives the sixth order polynomial as they got; on the other side, the same set of equations 
gives the fifth order polynomial as Kumar et al. (2006) obtained. The situation appears 
to be non-scientific, as the system of equations is a linear one. These are simple algebraic 
equations where the variables are to be eliminated. However, it is a matter of surprise that 
by solving these equations, two scientific groups are getting polynomials of different degrees. 
In the present discussion, we have attempted to short out this discrepancy. 

1 Introduction 

For application of magnetohydrodynamics (MHD) in solar physics as well as in plasma physics, 
dispersion relation, where uj is expressed as a function of k, plays a key role. A controversy for 
the degree of the polynomial in uj for dispersion relation appeared when Kumar et al. (2006, 
henceforth KKS) raised a point about the degree of the polynomial obtained by Porter et al. 
(1994, henceforth PKS) and by Dwivedi & Pandey (2003). Consequently, the results of PKS 
as well as of Dwivedi & Pandey (2003) were kept before a question mark. Energy equations of 
Dwivedi & Pandey (2003) was found erroneous (Klimchuk et al., 2004). After making correction 
in their energy equation, Dwivedi &: Pandey (2006, henceforth DP) made an attempt to show that 
the dispersion relation must be a sixth order polynomial. Since the results of an investigation 
involving MHD depend on the dispersion relation, it is important to resolve this controversy. 
This communication is an attempt to show that for the basic equations considered by DP, PKS 
and KKS, the dispersion relation comes out to be the same. In each case, it is a fifth order 
polynomial having the same coefficients. 

2 Basic equations of DP 

The basic equations used by DP as well as PKS are 




(1) 
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^ = Vx{vxB) (3) 

-^+1P{V. V) = {^-1) [Qth + Qvis - Qrad] (4) 

V = (5) 

rup 

with Qth = V. kVT and = ^ + .V). These equations are, respectively, the equation 
of continuity, equation of momentum, induction equation, energy equation and the equation of 
state. Here, p, v, ks, rUp, p, B, 7, T and 11 are, respectively, the total mass density, velocity, 
Boltzmann constant, proton mass, total pressure, magnetic field, ratio of the specific heats, 
temperature and the viscous stress tensor. 

For small perturbations from the equilibrium (PKS, KKS): 

P = Po + Pi 'v = vi B=Bo + Bi 

p = Po + Pi T = To + Ti n = no + ni 

where the equilibrium part is denoted by the subscript "0" and the perturbation part by the 
subscript "1". For the magnetic field taken along the z-axis, {i.e., Bo= Bqz) and the propagation 
vector k = kxX + kzZ, the equations (1) — (5) can be linearized in the following form. 

^ + Po(V.7i) = (6) 

Po^ = -Vpi + —(Vx5i)x 5o-V.no (7) 
at 47r 

^ = V X (^Ti X Bo) (8) 

^ + 7i?o(V. ^i) + (7 - = (9) 

^ = ^ + ^ (10) 
Po Po To 

For the perturbations that are proportional to exp[i{k ■ r —ojt)], equations (6) — (10) reduce to 
the following algebraic equations 

ujpi - po{kxVix + kzVi-^) = (11) 
ojpovix - kxPi - ^{kxBu - kzBix) + ^-^{klvi^ - '^hk^vu) = (12) 

i^Poviy + —{kzBiy) = (13) 

ivpoviz - kzpi + ^(4/c^t;u - 2kxkzVix) = (14) 

LdBix + k-^BoVix = (15) 

ujBiy + kzBoViy = (16) 

uBiz - kxBoVix = (17) 
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(18) 
(19) 



iojpi - ipoc^{k:rVix + kzVu) - (7 - l)'*!!^^^! 
PO Po To 

where = 7Po/Po- Equations (13) and (16) for the variables viy and Biy are decoupled from 
the rest and describe Alfven waves. The rest of the equations for pi, pi, Ti, Bi^, Bi^, vix and 
Viz describe damped magnetoacoustic waves. For elimination of the variables pi, pi, Ti, Bi^, 
Biz, vix and viz, we have 

:0 



(20) 






LaJ 













-Pokx 


-Pokz 


~kx 








j^kz 






(copo + '^kl) 




—kz 














2"?o I. I. 
3 ^z^^x 


{upo + ^kl 











LO 







kzBo 

















LO 




-kxBo 





iuj 





-il-^Hkl 










-ipoclkx 


-ipoclkz 


j_ 

PO 


_ J_ 

PO 


1 

To 

















3 Basic equations of KKS 



The basic equations used by KKS are 

dp 
dt 

Dv 



+ V.{pv) = 
1 



Dt 



Vp + i-(Vx B)y. B -v.n 
47r 



dB ^ ,^ 

_ = Vx (t; X S) 



Dp ^pDp 
'm ~ pDt 
2pkBT 



(7 — l)[Qt/i + Qvis — Qrad] 



p 



rUn 



(21) 
(22) 

(23) 
(24) 
(25) 



Comparison of the two sets [equations (1) — (5) and equations (21) — (25)] show that there is 
difference on the left side in the induction and energy equations. For this set of equations (21) 
— (25), after going through the same procedure as discussed in the preceding section, we get 
the equations (KKS) 

t^Pi - PoikxVix + kzViz) = (26) 

(27) 

(28) 





i^povix - kxPi - ^{kxBiz - kzBix) + ^-^(klvix - IkxkzVu) = 
47r o 



i^PoViy + — {kzBiy) 
47r 







i^PoViz - kzPi + —{4kzViz - 2kxkzVix) = 

LoBix + kzBoVix = 

ijBiy + kzBoviy = 

ujBiz - kxBoVix = 

iujpi — ipiujc^ — {1 — l)'^||^z^i = 



Pi 
Po 



Pi 
Po 



Ti 
To 



(29) 
(30) 
(31) 
(32) 
(33) 

(34) 
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Obviously, equation (18) is different from (33). DP started from the set of equations (1) — (5), 
but the results reported by them are as given here in equations (26) — (34). It appears that 
DP did not derive their equations (11) — (19), but adopted directly from KKS. It is noticeable 
that the present equations (11) — (19) are not available in the paper of PKS. Equations (28) 
and (31) for the variables viy and Biy are decoupled from the rest and describe Alfven waves. 
The rest of the equations for pi, pi, Ti, Bi^, Bi^. I'lx and vi^ describe damped magnetoacoustic 
waves. For elimination of the variables pi, pi, Ti, Bi^, Bi^, vix and viz, we have 
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4 Discussion and conclusion 

In order to resolve the controversy, now, we are left with two determinants (20) and (35), 
which correspond to the sets of equations used by DP (also PKS) and KKS, respectively. It is 
interesting to find out that both these determinants reduce to the following common determinant: 



kxiv u^po + i^kl-vlpok^ -^kxkz 

cqu - iij? copokx - ipQ(?gkxU copok^ - ipoc^kzUj 







where cq = (7 — 1)k|| fc^To/po and va = -Bo/v^Trpo • This determinant can be solved to get the 
following dispersion relation. 

Lo^ + iAuj"^ - Blo^ - iCu^ + Du + iE = Q 

where 
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3po 
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A = c^ + ^{kl + Akl) 
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+ v\cQk^ + 



Po Po 

U2i,2 A^.„..,2 U21.2 



A.ri{)v\k1k'^ 

3/30 



^ ^ ScQporjok^kj ^ Arjocovj^^kik' ^ v\c^k^k^ 



pI 



3/90 



2 u2 



E 



vj^copok^k 
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Hence, the dispersion relation obtained from both sets of the basic equations of DP (also PKS) 
and KKS is a fifth order polynomial in lo. The coefficients obtained here are the same as obtained 
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by KKS. It may finally be resolved that the dispersion relations derived by KKS for the basic 
set of equations (21) — (25) is correct. 

Though both the sets of basic equations produce a common dispersion relation, some points 
regarding the discrepancy between the induction and energy equations of the two sets can be 
noted as the following. 

In the induction equation (3) , the term DB /Dt can be expressed as 

DB dB ,^ 

Linearization of this equation gives 

DBx dB 



- + {vi .V) Bi 



Dt dt 

The second term on right side can be dropped as it is a product of two perturbations. Thus, we 
have 

DBi dBi 



Dt dt 

and the induction equation in the two sets give the same final equations. 

In the energy equation (4), the term 7p(V. v ) reduces to —ipocl{k,j.vix + kzViz) whereas in the 
equation (24), the term —('jp/p)(Dp/Dt) reduces to —ipiLoci.. However, after the calculations, 
no difference is found in the expression for dispersion relation. 
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